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Observation of hard scattering in 
photoproduction events with a large rapidity 

gap at HERA 



ZEUS COLLABORATION 



Abstract 

Events with a large rapidity gap and total transverse energy greater than 5 GeV have 
been observed in quasi-real photoproduction at HERA with the ZEUS detector. The dis- 
tribution of these events as a function of the 7p centre of mass energy is consistent with 
diffractive scattering. For total transverse energies above 12 GeV, the hadronic final states 
show predominantly a two-jet structure with each jet having a transverse energy greater 
than 4 GeV. For the two-jet events, little energy flow is found outside the jets. This 
observation is consistent with the hard scattering of a quasi-real photon with a colourless 
object in the proton. 
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1 Introduction 



In a recent publication [|l|, it has been shown that in photoproduction at HERA energies the 
contribution of diffractive processes account for about 36% of the total photoproduction cross 
section. Diffractive processes are generally believed to proceed via the t-channel exchange of 
a colour-singlet object, with vacuum quantum numbers and which carries energy-momentum, 
called the pomeron. The true nature of the pomeron is still far from clear. Ingelman and Schlein 
assumed that the pomeron emitted from the proton behaves like a hadron and suggested that 
it could have a partonic substructure which could be probed by a hard scattering process. The 
UA8 experiment at CERN later observed events containing two high-p^ jets in pp interactions 
tagged with leading protons (or antiprotons) 0]. These observations could be explained in 
terms of a partonic structure in the pomeron. 

In previous publications 0, |, ||, |^ evidence has been presented for events with a large 
rapidity gap in deep inelastic scattering (DIS). The energy dependence of the event rate as 
well as the approximate scaling present in the data pointed to a diffractive process of a leading 
twist nature. The hadronic final state of these events exhibited a small, but significant, rate 
for two-jet production in the 7*p frame. The characteristics of these events are consistent with 
an interaction between a virtual photon and partons in a colourless object from the proton. 
In this context the term 'pomeron exchange' is used as a generic name to describe the process 
which is responsible for creating events with a large rapidity gap. 

In this paper we report the observation of events with a large rapidity gap in a sample of 
events with high transverse energy produced in the photoproduction regime (Q^ ~ 0, where 
—Q^ is the four-momentum transfer squared carried by the virtual photon). The analysis of 
the associated hadronic final states in these events includes the study of jet structure to search 
for a hard scattering process in diffractive photoproduction. 

2 Experimental setup 

2.1 HERA machine conditions 

The data presented here were obtained with the ZEUS detector during the 1993 running period 
at the electron-proton collider HERA when 84 bunches of electrons with energy = 26.7 GeV 
collided with 84 bunches of protons of energy Ep= 820 GeV. In addition, 10 electron and 6 
proton non-colliding bunches were used for studies of beam induced background. The electron 
and proton beam currents were typically 10 mA. 

2.2 The ZEUS detector 

ZEUS is a multipurpose magnetic detector whose configuration has been described elsewhere 
1^. For the present analysis we used only some of the components within ZEUS. Charged 
particles are tracked by the vertex detector (VXD) |Q and the central tracking detector (CTD) 
which operate inside a thin superconducting solenoid providing an axial magnetic field of 1.43 
T . The solenoid is surrounded by a high resolution uranium-scintillator calorimeter divided into 
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three parts, forward (FCAL) covering the pseudorapidity Q region 4.3 > rj > 1.1, barrel (BCAL) 
covering the central region l.l > r] > —0.75 and rear (RCAL) covering the backward region 
—0.75 > rj > —3.8. The resulting solid angle coverage is 99.7% of An. The calorimeter parts 
are subdivided into towers which in turn are subdivided longitudinally into electromagnetic 
(EMC) and hadronic (HAC) sections. The sections are subdivided into cells, each of which is 



viewed by two photomultiplier tubes. The calorimeter is described in detail elsewhere |10]. For 
measuring the luminosity as well as for tagging the scattered electron in small processes. 



we use two lead-scintillator calorimeters ||TT]|. For these 'tagged' photoproduction events, the 
resulting values are less than 0.02 GeV^. 



2.3 Trigger conditions 

The data were collected with a three-level trigger. The First Level Trigger (FLT), based on 
a deadtime-free pipeline, selects inclusive photoproduction events with a calorimeter energy 
trigger. For FLT purposes, the calorimeter is subdivided into 896 trigger towers, each tower 
consisting of an EMC and a HAC segment. Events for this analysis were accepted by the FLT 
if the energy sum of EMC towers in the BCAL exceeded 3.4 GeV, or in the RCAL (excluding 
the cells adjacent to the beam pipe) exceeded 2.0 GeV, or the EMC towers exceeded 3.75 GeV 
in the RCAL towers adjacent to the beam-pipe. 

The Second Level Trigger (SLT) used information from a subset of detector components to 
differentiate physics events from backgrounds. The SLT rejected proton beam-gas background 
by timing measurements in the calorimeter cells; this algorithm reduced the RCAL and BCAL 
FLT rates by approximately 90% and 50%, respectively. 

The Third Level Trigger (TLT) had available the full event information on which to apply 
physics-based filters. The TLT applied stricter cuts on the event times and also rejected beam- 
halo muons and cosmic muons. The logic of the filter used in this analysis is described in the 
data selection section. 



3 Kinematics of photoproduction events 

In electron-proton scattering, photoproduction can be studied in the limit of small four-mo- 
mentum transfer, q, carried by the virtual photon, 7*. The kinematic variables used to describe 
the reaction e (k) + p (P) e {k') + X, are the following: the square of the total ep centre 
of mass energy: s = {k + P)^ ^ AEpE^ = 87576 GeV^; the four-momentum transfer squared 
carried by the virtual photon: Q"^ = —q^ = —{k — k')^; the Bjorken variable describing the 
energy transfer to the hadronic system: y = and the centre of mass energy squared of the 
7*p system: W'^ = {q + PY = ys. 

For events with a detected electron, the variable y can be obtained from 

Kl^^ (1) 



^The ZEUS coordinate system is defined as right-handed with the Z axis pointing in the (forward) proton 
direction and the X axis pointing horizontally towards the centre of the HERA rings. The pseudorapidity ry 
is defined as — ln(tan |), where the polar angle 9 is taken with respect to the proton beam direction from the 
nominal interaction point (IP). 
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where E'^ denotes the scattered electron energy and 9'^ the electron scattering angle. Alterna- 
tively, y can be determined approximately from the hadronic system using the Jacquet-Blondel 
expression 

yjB = — ^7-^^ — (2) 

with the sum running over all calorimeter cells i associated with the hadronic system. We denote 
by Ei the energy in cell i, pzi = Ei-cos9i and 9i is the angle of the centre of the cell with respect 
to the event vertex. The invariant mass, Mx, of the hadronic system detected in the ZEUS 
central detector is determined from the calorimeter cell information: Mj- = Ej^ — p\j, where 
Eh and ph are the energy and momentum of the hadronic system. Cells in the electromagnetic 
(hadronic) calorimeter sections with energies below 60 MeV (110 MeV) were excluded in the 
present analysis. 



4 Data selection 

The cuts to select hard photoproduction events are similar to those described in our previous 
publications 0, 0. In particular, we 



• reject events with an electron found in the calorimeter with an energy greater than 5 GeV 
and i/e < 0.7, to remove DIS background; 

• require 0.05 < yjB < 0.8, to remove beam-gas and those DIS events where the scattered 
electron was not identified and therefore was mistakenly included in the hadronic system; 

• reject events with an energy in the rear section of the calorimeter Ercal > 30 GeV, to 
remove background from non-ep collisions; 

• reject events with a missing transverse momentum, > 10 GeV/c; 

• require a vertex with —35 < Z <20 cm and a radial distance from the beam line R < 4 
cm; 

• require y .^ ^' < 0.9, at least two oppositely charged tracks with pr > 0.5 GeV, and 

one of the following: transverse energy in a cone outside 10° of the forward direction in 

excess of 5 GeV; or yjB > 0.28; or " < 0.8; or an electron with an energy larger 

than 5 GeV detected in the luminosity monitor; these criteria were essential to enable us 
to reach lower transverse energies than in previous publications JT^ ; 

• require a total transverse energy E^ > 5 GeV. 

The last condition in particular defines 'hard photoproduction' in the context of this anal- 
ysis. 

From an integrated luminosity of 0.55 pb~^, a sample of 417081 events passed these cuts. 
This sample has less than a 0.1% contamination from beam-gas interactions, as determined from 
the number of events originating from non-colliding electron and proton bunches. The cosmic 
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ray background, estimated from the rate of events outside of ep bunch crossings, is neghgible. 
The requirement that no electron is found in the ZEUS calorimeter ensures that < 4 GeV^. 
Monte Carlo studies, using the ALLM prescription which provides a smooth interpolation 
from deep inelastic scattering to = 0, show that for the accepted photoproduction events 
the median is 10~^ GeV^. The same Monte Carlo program predicts that 27% of the events 
should have a scattered electron measured in the electron calorimeter of the luminosity detector, 
in agreement with the observed fraction of 26.7%. 



5 The Monte-Carlo simulation 

The hadronic final states observed in hard photoproduction can be understood as the result of 
two different leading order (LO) mechanisms: 

• direct processes where the photon interacts with a gluon [g) in the proton giving rise to 
a quark-antiquark pair (Boson Gluon Fusion) or with a quark (g), generating a qg final 
state (QCD Compton), and 



resolved processes |JT5|| where a parton in the photon interacts with a parton in the proton. 



Calculations of the relative importance of these two competing mechanisms at HERA ener- 
gies have shown that the hard photoproduction cross section is dominated by resolved processes 
TE\. These expectations have been recently confirmed experimentally |T^, |TB|, |T^. 



The Monte Carlo generator PYTHIA 5.6 [^] was used to model standard hard photoproduc- 
tion processes. In this generator, the direct and resolved photon processes are each simulated 
using leading order matrix elements, with the inclusion of initial and final state parton showers. 
The lower cut-off on the transverse momentum of the generated final-state partons ptmin was 
chosen to be 2.5 GeV/c. The photon parton distributions were parametrized according to GRV- 
LO [0 while for the proton MRSD'_ was used. The Weizsacker- Williams approximation 



was used to describe the photon flux at the lepton vertex. This Monte Carlo simulation does 
not contain any explicit contribution from diffractive 7p interactions. 



Diffractive processes were simulated using POMPYT |21|. This is a Monte Carlo model 
in which, within the framework provided by PYTHIA, the proton emits a pomeron whose 
partonic constituents subsequently take part in a hard scattering process with the photon or 
its constituents. This model incorporates approximately energy independent cross sections as 
experimentally determined in hadron-hadron collisions. 

The photon contribution contains both direct and resolved processes and the parton density 



distribution is parameterised according to DC ||2^. The parton momentum densities of the 



pomeron are parameterised according to the hard distribution [Q, 21 

(3f{(3) = constant ■ P{1 - (3) (3) 

where (3 = Xparton/pom denotes the fraction of the pomeron momentum involved in the scat- 
tering.0 Two possibilities have been considered: one in which the partons in the pomeron 

^The soft parton density has not been considered since it does not describe our DIS data a. 
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are quarks (quarkonic pomeron) and a second one in which the partons are gluons (gluonic 
pomeron). The two predictions from the POMPYT model will be shown separately. 

In our previous publication on the study of jet production in DIS events with a large 
rapidity gap, the Nikolaev-Zakharov (N-Z) model |^ gave similar results to POMPYT. For the 



photoproduction processes studied here, the N-Z model again gives similar results to POMPYT. 

For all of the comparisons shown below, the number of Monte Carlo events has been nor- 
malised to the data in each figure and so only the shapes of the distributions may be compared. 
The Monte Carlo events were passed through reconstruction and selection procedures identical 
to those for the data. 



6 Results 



6.1 Events with large rapidity gaps 

Following [|, we define "^max as the maximum pseudorapidity of all calorimeter condensates in 
an event, where a condensate is defined as an isolated set of adjacent cells with summed energy 
above 400 MeV. The pseudorapidity of the condensate is then calculated from the angle of the 
energy weighted centre of the condensate with respect to the measured IP. The distribution of 
Vmax is shown in Fig la. The data presented in this and all following figures are not corrected for 
effects from detector acceptance and smearing. The dip in the rjmax distribution at rjmax ~ 1.1 
is a detector effect. Values of rjmax > 4.3, which are outside the acceptance of the calorimeter, 
occur when energy is deposited in many contiguous cells around the beam pipe in the forward 
(proton) direction. This region is sensitive to the fragmentation of the proton remnant which 
at HERA energies is not yet fully understood. The bulk of the events cluster around r]max ~ 4 
in fair agreement with the expectations of PYTHIA. In addition to this region of large rjmax a 
second class of events with ijmax < 1.5 is seen in the data. There are 6678 events with ijmax <l-5, 
corresponding to 1.6 % of the total sample. The beam-gas contamination in this subsample is 
less than 1% and that due to cosmic rays is 1.3%. 

The shape of the distribution for 'q^ax < 1.5 is not accounted for by standard Monte 
Carlo simulations for hard photoproduction processes as in PYTHIA. It is, however, in good 
agreement with the predictions of the POMPYT model, as illustrated in Fig. la. Note that 
the normalisations of the POMPYT and PYTHIA samples in this figure have been fixed to 
the number of data events below and above r]max = 1-5, respectively. According to POMPYT, 
for events with Et > 5 GeV the acceptance after the ijmax < 1-5 cut is about 9% (10%) for a 
quarkonic (gluonic) pomeron. 

In Regge phenomenology for soft processes the amplitudes for two-body scattering by 
pomeron exchange are characterized by a constant or slowly rising dependence on W, the 
centre of mass energy, while reggeon exchange leads to a power law decrease. Hence, the W 
dependence of the rate of large rapidity gap events is a sensitive measure of the type of exchange 
contributing to the scattering process. It is important to study whether the same phenomenol- 
ogy exists for hard scattering, the subject of this paper. The distributions inW = {yjB ■ s)^ for 
all events and for those with rjmax <1.5 are shown in Fig. lb. Also shown are the POMPYT 
predictions for events with r]max < 1-5. The good agreement of POMPYT with the data is 
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consistent with the assumption that the dominant mechanism for large rapidity gap events is 
pomeron exchange. 

The mass of the hadronic system was measured with the calorimeter as described at the 
end of section 3. The distribution of Mx for events with a large rapidity gap is shown in 
Fig. Ic. The resolution in the determination of Mx is between 10-20%. According to Monte 
Carlo calculations the reconstructed Mx values are typically underestimated by about 20% 
with respect to the generated masses. The Mx distribution is steeply falling for Mx values 
above 12 GeV. This is similar to the behaviour observed in DIS for events selected with similar 
cuts 1^. Neither the gluonic nor quarkonic POMPYT model gives a satisfactory description of 
the data. The correlation between Mx and r]max is displayed in Fig. Id. Note that, according 
to POMPYT, there are diffractive events at high Mx values which are suppressed by the rj^ax 
cut. 

If large rapidity gap events are interpreted as being due to pomeron exchange, then the 
fraction of the proton momentum carried by the pomeron Xpom/p can be determined from the 
mass of the hadronic system via the relation 

- 

Xpom/p — 'yy^- \^) 

For events with r^^ax < 1-5, Xpom/p clusters (not shown) around values of ~ 3 • 10~^ due to the 
applied cuts. The Xpom/p values are limited on the lower side to T-IO""^ by the Et > 5 GeV 
requirement while the "^max < 1-5 cut suppresses values of Xpom/p above 10~^. Similar Xpom/p 
values have been observed in DIS 0. 



6.2 Jet structure 



Evidence for multijet structure in hard photoproduction at HERA has been presented in |I2| , |T^, 
p!7|] . In the events with Et > 5 GeV a search was performed for jet structure using a cone-based 
jet algorithm in pseudorapidity (rj), azimuth (0) space subject to the Snowmass convention 
p5| . The cone radius R = (A^^ -|- A?]^) 2 in the algorithm was set to 1 unit. In order to ensure 
that for standard (non-diffractive) hard photoproduction the results are not biased by fragments 
from the proton remnant, whose fragmentation properties at these energies are not well known, 
calorimeter cells with polar angles smaller than 9° (r/ > 2.5) in the laboratory were excluded. 

In the first step of the jet algorithm, each calorimeter cell with a transverse energy in excess 
of 300 MeV is considered as a seed for the search. These seeds are combined if their distance 
in 77 — space, R, is smaller than 1 unit. Then a cone of radius i? = 1 is drawn around each 
seed and the calorimeter cells within that cone are combined to form a cluster. The axis of 
the cluster is defined according to the Snowmass convention: 7^'='"**^'' ^(pciuster^ ^j^g transverse 
energy weighted mean pseudorapidity (azimuth) of all the calorimeter cells belonging to that 
cluster. A new cone of radius 1 unit is then drawn around the axis of the cluster. All cells 
inside the cone are used to recalculate a new cluster axis. The procedure is iterated until the 
content of the cluster does not change. 

The energy sharing of overlapping clusters is then considered. Two clusters are merged if 
the common transverse energy exceeds 75% of the total transverse energy of the cluster with 
the lowest transverse enrgy; otherwise two different clusters are formed and the common cells 
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are assigned to the nearest cluster. Finally, a cluster is called a jet if E^^^^"^^ > 4 GeV. The 
transverse energy weighted mean pseudorapidity {rjjet) and azimuth {4>jet) were evaluated and 
jets were accepted for the present analysis if r/jet ^ 2, corresponding to polar angles larger than 
15°. 

We find that 91.4%, 6.5%, 2.0% and 0.1% of the events with total transverse energy > 
5 GeV belong to the zero-, one-, two- and three- or more-jet categories in the sample of events 
with rimax < 1-5. The fraction of zero-jet events depends strongly on the jet definition, in 
particular on the requirement that the jet transverse energy be larger than 4 GeV. 

To determine the background present in the two-jet sample with a large rapidity gap, we 
studied the empty as well as the non-colliding electron and proton bunches and set an upper 
limit of 1% for the beam-gas and cosmic ray backgrounds. For further confirmation that this 
sample of two-jet events is from photoproduction processes, we measured the fraction of two- 
jet events where the scattered electron is tagged in the electron calorimeter of the luminosity 
monitor with an energy between 5 and 22 GeV and for which < 0.02 GeV^. This was found 
to be 24% in agreement with both Monte Carlo expectations and the value for the complete 
high Et photoproduction sample. 

Figure 2a displays the rjmax distribution for all events with two or more jets. A sample 
of 132 two-jet events with rjmax < 1-5 is observed, corresponding to 0.63% of the hard photo- 
production two-jet sample. This number of events is not accounted for by the standard hard 
photoproduction processes as modelled by PYTHIA which predicts fewer than 18 ± 7 two-jet 
events with ijmax < 1-5. Therefore, the two-jet events observed with rj^ax < 1-5 are not just 



the tail of 'standard' hard photoproduction with two jets [12|. Both POMPYT samples give a 



good representation of the shape of the r]max distribution for these events. 

For events with 'qmax < 1-5, the distribution of the total transverse energy per event, Et, 
is shown in Fig. 2b for all events (open histogram), for events with at least one jet in the final 
state (cross hatched) and for those events with two or more jets (shaded). For E-r > 8 GeV, 
which is the minimum E^ for which two-jet production with E^ > 4 GeV is possible, 9.9% 
of the events are of the two-jet type. For Et > 12 GeV the majority of the events are of the 
two-jet type; an example of which is shown in Fig. 2c. For two-jet events the distribution of 
the transverse jet energies reaches up to values of 10 GeV, as shown in Fig. 3a. In Fig. 3b 
the difference in azimuth (A^-'^*) of the two jets in the large rapidity gap sample is displayed. 
The two jets are preferentially back-to-back in the transverse plane. The predictions from both 
POMPYT models are in fair agreement with these data. 



6.3 Transverse energy flow around the jet axis 

The profiles of the jets observed in high Et photoproduction events are compared for events 
with and without a large rapidity gap. Figures 3c, d show the transverse energy flows around 
the jet axis. To reduce the bias from the ij^ax < 1-5 cut, the jets for both samples were 
restricted to the region rjjet <0 and calorimeter cell energy deposits with rjceii > 1-5 were 
excluded. Fig 3c shows that the weighted azimuthal distributions (A0 = (j)ceii — 4>jet axis) 
in the jet core are essentially the same for the large rapidity gap events {rjmax < 1-5) and for 
the hard photoproduction sample with r]max > 1-5. However, the transverse energy outside of 
the jet core is about a factor of two larger for the r]jnax > 1-5 sample than for the large rapidity 
gap events. 
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Figure 3d shows the E'-r-weighted distribution for Arj, the difference in rapidity of a given 
calorimeter cell and the jet axis. For this figure, only energy deposits in the hemisphere con- 
taining the jet are included. The transverse energy flow in the core of the jet is the same for 
events with and without a large rapidity gap, indicating that the jets themselves are similar. 
However, jets produced in events without a large rapidity gap show at large (forward) A77 val- 
ues a significant flow. In contrast, the jet profile for the large rapidity gap sample is more 
symmetric and shows a factor of two to three less Ej- flow in the forward direction. These 
observations are similar to those presented in our previous publications on DIS events with and 
without a large rapidity gap [|, Q where it was demonstrated there is a suppression of colour 
flow between the outgoing nucleon system and the struck parton in events with a large rapidity 
gap. 



6.4 Momentum fractions 

We conclude that we have observed large rapidity gap events containing two jets, consistent 
with a photon pomeron hard scattering process leading to a two-jet final state. In 2 ^ 2 
parton scattering the momenta of the incoming partons can be calculated from the two partons 
in the final state. Let and Xp be the fractions of the photon and proton momenta carried by 
the initial state partons. If the final state partons are approximated by the two jets observed, 
energy-momentum conservation leads to the relationship: 

/_^yJ-J yzjhadrons 
_ nE+Pz)jets . . 

2Ep ■ 

In these expressions 'Jets' refers to a sum over all calorimeter cells that comprise the jets 
according to the cone algorithm employed. Direct photon processes are characterised by x^ 
values close to 1. 

Equations (5) and (6) were used to calculate x^ and Xp. The resolutions are comparable 
to those reported in an earlier publication on the hard photoproduction of two jets[|l^] and 
are in the range (10-20)%. The x^ distribution shown in Fig. 4a peaks near x^ = 1. The 
sum of the direct and resolved photon contributions as calculated by POMPYT for a gluonic 
pomeron gives a reasonable description of the data, as does a pure direct photon interaction 
with a quarkonic pomeron. In contrast, the shape of a purely resolved photon contribution is 
inadequate to describe the data, as seen by the dashed-dotted histogram in Fig. 4a which was 
calculated for a quarkonic pomeron with Eq. (3). 

The distribution of Xp, the fraction of the proton's momentum participating in the hard 
scattering, is displayed in Fig 4b and shows that most of the two-jet events have Xp values 
between 3-10^^ and 10~^. The low and high limits are due to the Et and rjmax cuts, respectively, 
as discussed in section 6.1. In this sample, the two jets populate the pseudorapidity region 
—2 < rjjet < 1 as shown in Fig. 4c, because the ijmax < 1.5 cut suppresses jets with tjjet > 1. 

In standard hard photoproduction at HERA, the centre-of-mass of the hard process is 
boosted in the proton direction and, therefore, in most two-jet events both jets go forward | |12[| . 
This is particularly true for the resolved photon process which is, in general, the dominant jet 
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production process at the transverse energies considered here. At centre-of-mass energies of 
W ~ 200 GeV, the direct photon process is expected to dominate only at transverse energies 



above 50 GeV [T^. In the selected sample of large rapidity gap events, Mx is an order of 
magnitude smaller than W, thus limiting the available phase space. The dominance of the 
direct component hence occurs at transverse energies an order of magnitude smaller than in 
standard hard photoproduction. Therefore the fact that two-jet events in the large rapidity gap 
sample are found to be mainly due to the direct coupling of the photon is well understood with 
the conventional parton density parametrisations of the photon. These features are reproduced 
by the POMPYT model. 



7 Conclusions 



We have observed photoproduction events with a large rapidity gap and large transverse energy 
at HERA. Their distribution as a function of the 7p centre of mass energy is consistent with a 
diffractive process. Hard scattering, with jets having transverse energies greater than 4 GeV, 
has been observed in these large rapidity gap events. For total transverse energies above 12 GeV 
the hadronic final state is dominated by two-jet production with the two jets being preferentially 
back-to-back in azimuth. For the two-jet events selected by the cuts used in this analysis, the 
fraction of the photon momentum participating in the hard scattering is close to one, suggesting 
that their production is dominated by direct photon processes. The two-jet events in the large 
rapidity gap sample show little energy outside the jet core. A natural interpretation of these 
events is the hard interaction of the photon with a colourless object inside the proton: the 
pomeron. These conclusions complement those drawn from the production of large rapidity 
gap events in deep inelastic scattering. 
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Figure 1: 

(a) The distribution of rjmax, the pseudorapidity of the most forward calorimeter condensate 
with energy above 400 McV for the photoproduction sample with Et > 5 GeV along with the 
predictions from PYTHIA (dashed line) and POMPYT with a quarkonic (solid line) or gluonic 
(dotted line) pomeron. 

(b) The distribution in W for all events and for those with rjmax < 1-5. The latter are well 
described by pomeron induced reactions as implemented in POMPYT. 

(c) The mass of the hadronic system Mx for events with a large rapidity gap as defined by 
Vmax < 1-5 along with the POMPYT predictions. 

(d) A scatter plot of the mass of the hadronic system, Mx, versus rjmax- 
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Figure 2: 

(a) The distribution of rjmax in the photoproduction sample with Et > 5 GeV and two or more 
jets along with the predictions from PYTHIA (dashed line) and POMPYT with a quarkonic 
(solid line) or gluonic (dotted line) pomeron. 

(b) The distribution of the total transverse energy Et for the photoproduction sample with a 
large rapidity gap and, in addition, for the subsample of those events with at least one (cross 
hatched ) and at least two (shaded ) jets in the final state. 

(c) A display of a two-jet event in hard photoproduction with a large rapidity gap. 
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Figure 3: 

(a) The transverse energy distribution of the jets in the two-jet sample in photoproduction with 
Vmax < 1-5 along with the POMPYT predictions for a quarkonic (solid hne) or gluonic (dotted 
line) pomeron. 

(b) The A0-^^* distribution for two-jet events in photoproduction with rjmax < 1-5 along with a 
comparison with the POMPYT predictions as in (a). 

(c, d) The transverse energy weighted profiles for jets with rjjet < in photoproduction events 
with Et > 5 GeV together with POMPYT predictions as in (a). Two samples are presented: 
one with rjmax <l-5 and the other with rj^ax > 1-5; (c) shows the azimuthal profiles and (d) 
shows the pseudorapidity profiles. 
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Figure 4: 

(a) The a;.y distribution for photoproduction events with a large rapidity gap (dots) compared to 
POMPYT predictions. The soUd hne represents the direct photon interaction with a quarkonic 
pomeron. The dotted hne denotes the sum of the direct and resolved interactions of the photon 
with a gluonic pomeron. Also shown are the results from POMPYT with a quarkonic pomeron 
using only the resolved photon processes (dash-dotted histogram). Each histogram has been 
normalised separately to the number of data events. 

(b) The measured fraction of the proton momentum, Xp, entering the hard scattering for two-jet 
events, along with the POMPYT predictions for a quarkonic (solid line) or gluonic (dotted line) 
pomeron. 

(c) The jet pseudorapidity distribution for the two-jet sample and for the POMPYT predictions, 
as in (b). 



